Photonuclear reactions of actinide and pre-actinide nuclei at intermediate energies 
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Photonuclear reaction is described with an approach based on the quasideuteron nuclear pho- 
toabsorption model followed by the process of competition between light particle evaporation and 
fission for the excited nucleus. Thus fission process is considered as a decay mode. The evaporation- 
fission process of the compound nucleus is simulated in a Monte-Carlo framework. Photofission 
reaction cross sections are analysed in a systematic manner in the energy range ~ 50-70 MeV for 
the actinides ^^^Th, ^^^U, ^^'^U, ^^*U and ^^^Np and the pre-actinide nuclei ^°**Pb and ^°^Bi. The 
study reproduces satisfactorily well the available experimental data of photofission cross sections at 
energies ~ 50-70 MeV and the increasing trend of nuclear fissility with the fissility parameter jA 
for the actinides and pre-actinides at intermediate energies [~ 20-140 MeV]. 
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I. INTRODUCTION 



In recent years the study of photofission has attracted 
considerable interest. When a gamma (photon) above 
the nuclear binding energy of an element is incident on 
that element, fission may occur along with the emission 
of neutrons and the number of neutrons released by each 
fission event is dependent on the element. With high 
enough photon energy it is possible to induce fission in 
most elements. For lighter elements the cross sections are 
quite small. There is also the difficulty in distinguishing 
fission from various fragmentation processes, particularly 
for light elements. These reasons practically limit all in- 
vestigations to the upper half of the Periodic Table. Re- 
cently the experimental techniques are improved by using 
the solid state track and parallel plate avalanche detec- 
tors which make it possible to study fission with very 
low yields. Moreover, the production of highly neutron- 
rich radioactive nuclei through photofission will open new 
perspectives to explore very exotic nuclei far away from 
the valley of stability, especially in the vicinity of ^^Ni. 

The use of the energetic electrons is a promising mode 
to get intense neutron-rich ion beams pj]. The reason is 
that although the fission cross section at the giant dipo- 
lar resonance energy is « 0.16 barn for ^38^ (against 
« 1.6 barn for 40 MeV neutrons induced fission) but the 
electrons/7-photons conversion efficiency is much more 
significant than that of the deuterons/neutrons. The en- 
ergetic beam (~ 50 MeV ) of incident electrons can be 
slowed down in a tungsten (W) converter or directly in 
the target (U), generating Bremsstrahlung 7-rays which 
may induce fission. In the photofission method for the 
(neutron-rich) radioactive ion beam production, nuclei 
are excited by photons covering the peak of the giant 
dipolar resonance. 
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Isotopes of plutonium and other actinides tend to be 
long-lived with half-lives of many thousands of years, 
whereas radioactive fission products tend to be shorter- 
lived (most with half- lives of 30 years or less). Many of 
the actinides are very radiotoxic because they also have 
long biological half-lives and are a emitters as well. In 
transmutation the intention is to convert the actinides 
into fission products. The fission products are very ra- 
dioactive, but the majority of the activity will decay 
away within a short time. From a waste management 
viewpoint, transmutation of actinides eliminates a very 
long-term radioactive hazard and replaces it with a much 
shorter-term one. Accelerated radioactive decay has been 
proposed by bombarding spent fuel with electromagnetic 
(photon) rays. The 7-rays have also been suggested to 
be very effective, but 7-rays are very difficult to produce 
and may need to be precisely tuned to the target actinide 
or fission product. 

The aim of the present work is to obtain photofission 
cross section. From the photonuclear reactions, it is pos- 
sible to gain insight into the fission mechanism, especially 
about how the barrier varies. Information can also be 
obtained about the excitation of a nucleus by absorp- 
tion of a high-energy quantum and the subsequent de- 
excitation process. At intermediate energies [~ 20-140 
MeV] such a reaction can be described as a two step pro- 
cess. The incoming photon is assumed to be absorbed by 
a neutron-proton [n-p] pair inside nucleus, followed by a 
mechanism of evaporation-fission competition. Photofis- 
sion cross sections of various nuclei are calculated using a 
Monte-Carlo method for the evaporation-fission calcula- 
tion and compared with the available experimental data 
at energies ~ 50-70 MeV. 



II. THE PHOTOABSORPTION MECHANISM 

The dominant mechanism for nuclear photoabsorp- 
tion at intermediate energies is described by the 
quasideuteron model [3] which is employed to evaluate 
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the total photoabsorption cross section in heavy nuclei. 
It is based on the assumption that the incident photon 
is absorbed by a correlated n-p pair inside the nucleus, 
leaving the remaining nucleons as spectators. Such an 
assumption is enforced when wavelength of the incident 
photon is small compared to the nuclear dimensions. The 
total nuclear photoabsorption cross section crj is then 
proportional to the available number of n-p pairs inside 
nucleus and also to the free deuteron photodisintegration 
cross section Ud{E^) and is given by : 



NZad{E^)fp{E^ 



(1) 



where is the incident photon enegy, N, Z and A 
are the neutron, proton and mass numbers respectively, 
L/A factor represents the fraction of correlated n-p pairs. 
Thus, Levinger's constant L is a factor which measures 
the relative probability of two nucleons being near each 
other in a complex nucleus compared with that in a free 
deuteron. The function fp{E^) accounts for the reduc- 
tion of the n-p phase space due to the Pauli exclusion 
principle. A systematic study of the total nuclear pho- 
toabsorption cross section data Q in the intermediate 
energy range shows that 

fp{E^) ^e-^/^^ where D = 0.72A"-^^MeY . (2) 

For photon energies upto the pion threshold, the Eq. (2) 
along with the damping parameter D provided above 
agrees reasonably well with the approach based upon 
phase space considerations Q using Fermi gas state den- 
sities that conserve linear momentum for the Pauli block- 
ing effects in the quasideuteron regime of hard photon 
absorption. 

The free deuteron photodisintegration cross section [5] 
is given by 



(Jd{E^) = 



61.2 [E-, - B]3/2 



E-^ 



lb] 



(3) 



where B=2.224 MeV is the binding energy of the 
deuteron. The quasideuteron model of nuclear photoab- 
sorption is used together with the modern root-mean- 
square radius data to obtain Levinger's constant 



i = 6.8- 11.2^ 



-2/3 



5.7A 



-4/3 



(4) 



of nuclei throughout the Periodic Table and is in good 
agreement @ with those obtained from the experimen- 
tally measured values. At the quasideuteron energy 
range, as a consequence of the primary photointerac- 
tion, 7-|-(n-|-p)^n*-|-p*, in most of the cases excited com- 
pound nuclei are formed with the same composition as 
target nucleus where both neutron and proton are re- 
tained inside the nucleus and the probabilities that either 
neutron escapes or proton escapes or both neutron and 
proton escape from the nucleus are extremely low. 



III. THE NUCLEAR EXCITATION AND 
FISSION PROBABILITY 

After absorption of a photon (7) with incident energy 
Ej (as measured in the laboratory frame), the nucleus 
with rest mass mo recoils with a velocity Vr given by 



E^c 



iriQC 



(5) 



which is also the velocity of the centre of mass Vcm of the 
7-nucleus system where c is the speed of light in vacuum. 
In the centre of mass frame, momenta of the 7 and the 
nucleus are the same and equal to Pcm- 



Pcm 



E-yiriQC 



(6) 



and the kinetic energy Er of the recoiling nucleus in the 
laboratory frame is 



Er 



J „2 



= E~^ + niQC — [moc {2E-y + toqC )] 



2^1/2 



(7) 



where toq and TOq are the rest masses of the nucleus be- 
fore and after the photon absorption respectively. 

The recoiling nucleus can be viewed as a compound nu- 
cleus having the same composition as the target nucleus 
but with the excitation energy 



moc 



moc^[(l 



1] 
(8) 



and in this case E* is also equal to the kinetic energy 
in the centre of mass frame (which is sum of the kinetic 
energies of the 7 and the nucleus moving in the centre 
of mass frame). The kinetic energy Er of the recoiling 
nucleus in the laboratory frame can now be rewritten as 
the obvious result 



Er — E^ — E"^ 



(9) 



This excited compound nucleus then undergoes succes- 
sive evaporation of neutrons and other light particles or 
fission. Thus, the fission is considered as a decay mode. 
The photofission cross section ct/ is a product of the nu- 
clear photoabsorption cross section crj and the total fis- 
sion probability (fissility) / and is, therefore, given by 



(10) 



The statistical approach for nucleon and light-particle 
evaporation and nuclear fission is an appropriate scheme 
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for calculation of the relative probabilities of different 
decay modes of the compound nucleus. Such statistical 
decay of the compound nucleus is the slow stage of the 
photonuclear reaction. According to the standard Weis- 
skopf evaporation scheme 7] , the partial width Tj for the 
evaporation of a particle j — n, p, ^H, '^H, '^He or ^He is 
given by 



E'-Bi 



<n.iE)p,{E*^B,~E)EdE 
(11) 

where Sj, /ij, Vj and Bj are the spin, reduced mass, 
Coulomb barrier and the binding energy of the parti- 
cle j, respectively. cr^„„(i?) is the cross section for the 
inverse reaction which means the capture reaction cross 
section of the particle j to create the compound nucleus. 
PCN and pj are the nuclear level densities for the ini- 
tial and final (after the emission of the particle j) nuclei, 
respectively. 

The Bohr- Wheeler statistical approach [1] is used to 
calculate the fission width of the excited compound nu- 
cleus. This width is proportional to the nuclear level 
density Pf at the fission saddle point: 



r 



2npcNiE*) Jo 



E*-Bf 



Pf{E* - Bf - E)dE (12) 



where Bf is the fission barrier height. The diffuse surface 
nucleus Sierk's [§] fission barriers (Bj) are used for these 
calculations. The decay of the excited compound nucleus 
[lO'l is simulated using the Monte-Carlo method [ll| . The 
competition between the various decay channels at each 
step of the evaporation chain is determined by the rela- 
tion between the partial widths for particle evaporation 
and fission, Eqs. (11) and (12), respectively. Finally, in 
order to calculate the fission probability, the total num- 
ber of fission events in a computer run is counted and 
divided by the total number of photoabsorption events. 
Evaporation from excited fission fragments is also taken 
into account. 



IV. COMPETITION BETWEEN LIGHT 
PARTICLE EVAPORATION AND FISSION 

The quantitative description of the process is based 
on the liquid drop model (LDM) for nuclear fission by 
Bohr and Wheeler P and the statistical model of nuclear 
evaporation developed by Weisskopf Q . Accordingly, the 
probability of fission relative to neutron emission is calcu- 
lated using Vandenbosch-Huizenga's equation [l3| given 

by 



Tf _ Koan[2{afE})^ -1] 



r„ 



iA-safE* 



exp[2[(a/i?;)^-(a„i?:)^]] 
(13) 



where E^ = E* - B„ and EJ = E* - Bf are the nu- 
clear excitation energies after the emission of a neutron 
and after fission, resepctively, where Bn is the binding 
energy of the emitted neutron. r„ and F / are the partial 
widths for the decay of the excited compound nucleus 
via neutron emission and fission, resepctively, and the 
parameters a„ and af are the level density parameters 
for the neutron emission and the fission, respectively and 
Kq = It? /2mrQ where m and tq are the neutron mass and 
radius parameter respectively. The emission probability 
of particle k relative to neutron emission is calculated 
according to the Weisskopfs statistical model 



r„ 



exp[2[(a.i?^)i-(a„i?:)i]] 



(14) 

where Uk is the level density parameter for the emis- 
sion of the particle fc, ^kHn — 1 for k — 2 for 
k =4He, 1 for k =2H, 3 for k and 2 for /c = ^He. 
Ek = E* — {Bk + Vk) is the nuclear excitation energy 
after the emission of particle k [3] ■ -Sfc are the binding 
energies of the emitted particles and Vk are the Coulomb 
potentials. The evaporation- fission competition of the 
compound nucleus is then described in a Monte-Carlo 
framework. Any particular reaction channel r is then 
defined as the formation of the compound nucleus via 
photoabsorption and its decay via particle emission or 
fission. Thus, fission is considered as a decay mode. The 
photonuclear reaction cross sections ar are calculated us- 
ing the equation ar = a'^rir/N where is the number 
of events in a particular reaction channel r and A'' is to- 
tal number of events that is the number of the incident 
photons. 



V. COMPARISON OF THE PHOTOFISSION 
CROSS SECTION ESTIMATES WITH 
MEASURED DATA 

Each calculation is performed with 40000 events using 
a Monte-Carlo technique for the evaporation-fission cal- 
culation. This provides a reasonably good computational 
statistics. The photofission cross sections are calculated 
at different energies for various elements. Results of these 
calculations corresponding to the available experimental 
data at 52 MeV, 69 MeV, 120 MeV and 145 MeV are 
listed in Table-F The statistical error in the theoretical 
estimates for the photofission cross sections are calcu- 
lated using the equation cr/ ± Act/ — a"[[nf ± yJrrf]/N 

which implies that Act/ — ^Ja^/jf/N . For the cases, 
where not a single fission event occured in N events, the 
upper limits of the cross sections are calculated using the 
equation ct/ = cr^ /N where N[= 40000] is the number of 
incident photons. 

The results of the present calculations should be com- 
pared with measured photofission cross sections with cer- 
tain degree of caution. In this regard, it is pertinent to 
mention here that the experimental data listed in Table-I 
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TABLE L Comparison between the measured and the calcu 
lated photofission reaction cross sections. 



E =52 MeV 



E =69 MeV 



Target 


7 


Expt. 


Calc. 


Calc. 


nuclei 




J 




(Tf 

J 




MeV 


mb [Ref.] 


mb 


mb 




52 


20±2 [14] 


17.6 


17.61± 0.09 




69 


19±4 [14] 


15.6 


15.59 ±0.05 


235y 


52 


16±4 [15] 


17.6 


17.58± 0.09 




69 


15±3 [15] 


15.5 


15.54 ±0.08 


238 -(J 


52 


14±2 [14] 


17.5 


17.46± 0.09 






16±2 [15] 










32±2 [16] 










25±3 [17] 








69 


15±1 [15] 


15.5 


15.52± 0.08 






23±1 [15] 










13±2 [14] 








52 


8.6±0.6 [15] 


17.5 


12.41 ±0.07 




69 


9±1 [15] 


15.4 


13.82 ±0.07 


209gj 


52 


(16±l)xl0~^ [16] 


17.4 


0.15 ±0.01 






(70±12)xl0~^ [18] 










(24±3)xl0-^ [19] 








69 


(8.0±0.6)xl0"^ [16] 


15.0 


0.26 ±0.01 






(18±3)xl0"^ [18] 








120 


0.20±0.06 [20] 


9.3 


0.87 ±0.01 




145 


0.31±0.08 [20] 


7.6 


1.21 ± 0.02 




52 


(1.9±0.3)xl0"^ [16] 


17.3 


0.039 ±0.004 






3x10"^ [21] 








69 


(12±2)xl0"^ [16] 


14.9 


0.097±0.006 






(18±3)xl0"^ [21] 






i^^Au 


120 


(59±38)xl0"^ [20] 


9.0 


0.073 ±0.004 




145 


(ll±4)xl0-^ [20] 


7.3 


0.117±0.005 


i«iTa 


120 


(13±4)xl0"^ [20] 


8.5 


0.0023±0.0007 




145 


(9.7±3.0)xl0-3 [20] 


6.9 


0.0035±0.0008 


174Yb 


52 


(3.2±0.5)xlO-'5 [16] 


16.7 


< 4.2x10"" 




69 


(6±l)xl0"'^ [16] 


13.9 


< 3.5x10"'' 


is^Sm 


69 


(1.8±4)xl0"^ [16] 


13.0 


< 3.3x10"'' 


51y 


120 


(76±25)xl0"^ [20] 


3.0 


0.0019±0.0004 




145 


(78±29)xl0~^ [20] 


2.3 


0.0033 ±0.0004 



10.0000 



for = 69 MeV, 120 MeV and 145 MeV correspond to 
the effective mean energies of 69 MeV, 120 MeV and 145 
MeV respectively, for the incident quasi-monochromatic 
beams of photons. However, those for = 52 MeV 
actually correspond to the incident beam energy of 52 
MeV monochromatic photons. The same holds for the 
cross sections for the evaporation residues calculated for 
these photon induced reactions. But the theoretical cal- 
culations are performed exactly at photon energies of 
52 MeV, 69 MeV, 120 MeV and 145 MeV. Neverthe- 
less, present calculations provide good estimates of the 
photofission cross sections for the actinides. In the pre- 
actinide to medium mass region the agreement with the 
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FIG. 1: The plots of cross sections aER as a function of mass 
number Ar of the evaporation residues for the pre-actinides 
2°**Pb and ^''''Bi at Ej = 52 MeV, 69 MeV. The fission cross 
sections of the fissioning nuclei are also shown (marked /) . 



measured photofission cross sections is reasonable. 



VI. RESULTS AND DISCUSSION 

The photonuclear reaction cross section calulations at 
intermediate energies within a Monte-Carlo framework 
for simulation of the evaporation-fission competition are 
performed assuming 40000 incident photons for each cal- 
culation which provide a reasonably good computational 
statistics. The cross sections of fission and evapora- 
tion residues (cr^fl) as a function of mass number Aji 
of the evaporation residues for the pre-actinides ^ospj^ 
and ^"^Bi at the incident photon energies of 52 MeV and 
69 MeV for the present calculations are shown in Fig.l. 
This is to show the relative probabilities of other evapo- 
ration processes compared to fission. The results of the 
fission cross sections and the fissility at Ej = 20, 40, 60, 
80, 100, 120 and 140 MeV using 40000 events for the 
Monte-Carlo calculations are tabulated in Table-II. The 
increasing trend of the nuclear fissility with the fissility 
parameter Z"^ /A for the actinides at intermediate ener- 
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TABLE II: Variation of the calculated photofission reaction cross sections of actinides with the incident photon energy (Ej) 
and the fissility parameter /A. 



Target /A Calculated ^^[MeV] E 


J [MeV] E 


J [MeV] E 


J [MeV] E 


T [MeV] E 


T [MeV] E 


J [MeV] 


nil ^^Idi 




fill "aTrfit^r 
Vc^ LldliLl L V 




40 


60 


80 


100 


120 


140 


239 1-1 


36.97 


(T„ (mb) 


10.19 


10.04 


16.78 


14.24 


11.97 


10.13 


8.68 






CT/ (mb) 


10.19 


10.04 


17.78 


14.24 


11.97 


10.13 


8.68 






f (%) 


99.98 


100.0 


100.0 


100.0 


100.0 


100.0 


99.99 


237at 

"^'Np 


36.49 


aa (mb) 


10.31 


18.06 


16.74 


14.18 


11.90 


10.07 


8.62 






CT/ (mb) 


10.24 


18.05 


16.73 


14.18 


11.90 


10.07 


8.62 






J" /Cr/ \ 

f (%) 


99.38 


99.99 


99.99 


99.99 


100.0 


99.99 


99.99 


233 -y 


36.33 


(T„ (mb) 


10.58 


10.15 


16.71 


14.11 


11.82 


9.99 


8.54 






(7/ (mb) 


10.40 


10.14 


16.70 


14.11 


11.81 


9.98 


8.53 






J" /Cr/ \ 

f (%) 


98.34 


99.95 


99.97 


99.99 


99.97 


99.96 


99.93 


235-Q 


36.02 


cr„ (mb) 


10.42 


18.07 


16.69 


14.12 


11.84 


10.01 


8.57 






cr/ (,mDj 




18.05 


16.68 


14.11 


11.84 


10.00 


8.56 






/ (%) 


95.83 


99.91 


99.97 


99.96 


99.99 


99.97 


99.95 


238 -[J 


35.56 


o-I (mb) 


10.19 


17.94 


16.66 


14.13 


11.87 


10.04 


8.60 






0-/ (mb) 


8.94 


17.76 


16.64 


14.12 


11.86 


10.04 


8.59 






/(%) 


87.71 


98.97 


99.87 


99.96 


99.95 


99.96 


99.97 


232 rpj^ 


34.91 


erf (mb) 


10.57 


18.05 


16.59 


13.99 


11.72 


9.89 


8.46 






0-/ (mb) 


2.20 


10.18 


13.86 


13.27 


11.45 


9.70 


8.27 






/ (%) 


20.87 


56.38 


83.53 


94.79 


97.76 


98.06 


97.83 



gies [~ 20-140 MeV] are observed. 

Present calculations provide excellent estimates of the 
photofission cross sections for the actinides except for 
^^^Th where it somewhat overestimates the fission cross 
sections. However, the general increasing trend of nu- 
clear fissility with fissility parameter is retained. For the 
pre-actinides, ^"^Pb and ^"^Bi, the photonuclear reaction 
cross sections show peaks for evaporation residues around 
Ah = 203 (Pb), 201 (Pb) and 204 (Bi), 203 (Bi) respec- 
tively, for the incident photon energies = 52 MeV, 69 
MeV. The total number of events (that is the number of 
incident photons) for each run limits the calculations of 
too low fission cross-sections. The fissility for thorium 
[2^ 1 and several uranium isotopes Q was found to be 
lower than that for neptunium, showing that nuclear fis- 
sility does not saturate for those nuclei, remaining at a 
value below hundred percent even at high incident pho- 
ton energies. The present theoretical study corroborates 
this behaviour. We have included ^^^Pu in our study, 
though it is not a naturally occuring element but is read- 
ily produced as spent fuel. 



VII. SUMMARY AND CONCLUSION 



In summary, the cross sections for the fission and the 
evaporation residues are calculated for photon induced 
reactions at intermediate energies. Monte-Carlo calcu- 
lations for the evaporation-fission competition are per- 
formed assuming 40000 incident photons for each calcu- 
lation. These many number of events for each calculation 
provide a reasonably good statistics for computationally 
stable results. Present calculations provide excellent esti- 
mates of the photofission cross sections for the actinides. 
In the pre-actinide to medium mass region the agreement 
with measured photofission cross sections is reasonable. 
The process of photofission of heavy nuclei is considered 
in terms of production of fission fragments and is pro- 
jected as a viable method for the production of neutron 
rich nuclei for radioactive ion beam (RIB). 
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